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Abstract 

We perform a NNLO (and a LO) global parton analysis in which we include the new 
precise data for deep inelastic scattering from HERA and for inclusive jet production at 
the Tevatron, together with the improved knowledge of the three-loop splitting functions. 
The results are compared with our recent NLO analyses. The LO fit produces significantly 
worse results in general, but gives a surprisingly good fit to the Tevatron high-E'i^ jet data. 
For the approximate NNLO analysis we notice a slight improvement in the quality of the 
global fit, and find that the partons are changed by up to 10% at = 10 GeV^, in 
particular in the x < 0.01 regime. 



t Royal Society University Research Fellow 



The global parton analyses of deep inelastic scattering (DIS) and related hard scattering data 
are generally performed at NLO order. Recently much effort has been invested in computing 
NNLO QCD corrections to a wide variety of partonic processes and therefore we need to 
generate parton distributions also at NNLO so that the theory can be applied in a consistent 
manner. Analyzing DIS at NNLO is important in itself in that we may be able to investigate 
the hierarchy LO NLO — > NNLO in the process where there are perhaps the most precise 
data. Recently, there has been a significant increase in the precision of the structure function 
data from the HERA experiments 0-0 and these, together with the inclusive jet production 
measurements at the Tevatron 0, |^, were used to generate an updated version of the NLO 
partons ip.Q 

To carry out the analogous analysis at NNLO we need both the relevant splitting functions 
as well as the coefficient functions. However, although the deep inelastic coefficient functions 
are known at NNLO the splitting functions are not yet available at this order. On the 



other hand, valuable partial information is known in the form of the few lowest moments of the 



splitting functions [12|, which tightly constrain the large x behaviour. When this is combined 



with the known most singular log 1/x behaviour [|T^ , it greatly limits the possible behaviour of 
the functions down to quite small values of x. Indeed, van Neerven and Vogt |jT^ constructed 
a range of compact analytic expressions for the splitting functions that are all compatible with 
this information. In particular, they provided closed expressions which represent the fastest 
and the slowest permissible evolution. In an exploratory NNLO global analysis of the data |TB 



we, in turn, used these expressions to obtain bands representing the spread of allowed parton 
distributions. For example, the range of uncertainty of the NNLO gluon distribution was found 
to be about ±30% (±15%) at x ~ 10-^(10-=^). 

Subsequently, the moment calculation of the three-loop splitting functions was extended to 
give two higher moments [|TB[. This has allowed van Neerven and Vogt [|T^ to provide analytic 
expressions for the NNLO splitting functions which considerably increase the reliability of the 
evolution in the region x ^ 10"'^. In other words, the allowed bands of permissible splitting 
functions has shrunk considerably, even, somewhat surprisingly, at small x^ The splitting 
functions are located near the 'slow evolution' edge of their previous bands. 

Here, we use these improved constraints on the splitting functions to perform a new NNLO 
global parton analysis of the available data sets, which include the new more precise structure 
function measurements from HERA and the inclusive jet data from the Tevatron ^. 

The analysis uses the same form of the starting parameterizations and the same data sets 
[H-H, pl-H, with the same and cuts, as our recent NLO global analysis In 



^As in 1^, we use the DO data obtained using the cone algorithm rather than those obtained jl^ using the 
kx algorithm, sinee the former have greater rapidity coverage and have been studied more extensively. 

^The allowed band contains the vast majority of permissible behaviours, but it is still possible to find perfectly 
acceptable forms which lie considerably outside this region, particularly at very small x] i.e. the bands should 
not be regarded as anything like one-sigma errors in a Gaussian distribution. 
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particular, we fit to DIS data with > 2 GeV^ and > 12.5 GeV^^ In going to a NNLO 
analysis, we use the three-loop expression for as, in the MS scheme, and we extend the matching 
conditions when evolving through the heavy quark thresholds. The details of the heavy quark 
prescription are very similar to those in [|T^ , but we add a brief discussion of important points in 
the Appendix. Although the deep inelastic data are described at NNLO, some other processes 
are not yet calculated to this accuracy, and so we have to use NLO expressions to fit to the 
data for (i) Tevatron jet production, (ii) the charged lepton rapidity asymmetry from W boson 
hadroproduction and (iii) the Drell-Yan process^. We note that the effect of the threshold 



corrections at NNLO to Tevatron jet production has been calculated in [|30[, and this reduces 
scale uncertainty and perhaps is a good approximation to the full NNLO correction. The result 
is approximately a 3 — 4% increase in normalization of the theory prediction (not surprising 
since the NLO correction is itself about a 7 — 10% raise, largely independent of Et), which 
would be completely subsumed in the normalization error of the data. 

The optimum global NNLO fit is obtained with the starting parameterizations of the partons 
at Ql = 1 GeV^ given by 

xuv = 0.262a;°=^^(l-x)^-^°(l + 3.83x°-^ + 37.65x) (1) 

xdv = 0.061x°-^^(l -x)^-°^(l + 49.05a;°-^ + 8.65x) (2) 

xS = 0.759x-°-i^(l-x)^-^'^(l-1.34a;°-^ + 7.40x) (3) 

xg = 0.669x°-°°(l-x)3-9^(l + 6.98x°-^- 3.63x) - 0.23a;-°-'^(l-a;)^-^ (4) 

The fiavour structure of the light quark sea is taken to be 

2u,2d,2s = 0.4^ -A, 0.4^ + A, 0.2^ (5) 



with s = s, as implied by the NuTeV data |T9[, and where 



xA = x{d-u) = 1.432x^-2^(1 -x)^-^^(l + 9.86x - 29.04x2). (6) 

The masses of the quarks are taken to be rric = 1.43 GeV and = 4.3 GeV, the former 
giving the best fit to the charm structure function data. The value of a5(M|) = 0.1155, and 
is displayed in Table 1. This fit corresponds to the average of two extreme choices for the sets 
of splitting functions in Ref. and we can explore the uncertainty of the NNLO evolution 
by instead choosing either extreme. These 'slow' and 'fast' evolutions we denote by A and B 
respectively. 

For the sake of studying the progression from LO through NLO to NNLO, we also perform 
a pure LO fit, for which the starting parameterization is 

xuv = 0A74:X°-^%l-xf^^{l-l.32x°-^^ + I9.56x) (7) 



•^As for the NLO fit, we have investigated the consequences of changing these cuts, and introducing cuts at 
small X. The results are similar to NLO, and will be reported in p9| |. 

■^For Drell-Yan scattering the inclusive rate is known at NNLO, but not the differential distribution in xp 
that we use in the analysis. 
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xdv = 0.668x°-^3(l-x)^-°3(l-0.83x°-^ + 7.68x) (8) 

xS = 0.458x~°-^^(l-x)^-^^(l + 0.025x°-^ + 7.63x) (9) 

xg = 3.08x°-^°(l -2.96x°-^ + 9.26x) (10) 

with 

xA = x{d-u) = 4.163x^-^^(1 - a;)^-^^(l + 7.20x - 24.8x2). (11) 

In this case asiM"^) = 0.130 (again see Table 1), and the parameterization of the gluon is 
simpler because the fit no longer requires a negative input gluon at small x at this order. 



Table 1: The QCD coupling, and corresponding A parameter for rij = 4 (in the MS scheme 
with = in all cases), for the LO, NLO and NNLO fits. Note that as we progress from 
order to order a very different value of Aqcd is needed to give the same as(M^), that is as{Q'^) 
is defined differently at different orders. 





as{Ml) 


Aqcd (MeV) 


LO 


0.130 


220 


NLO 


0.119 


323 


NNLO 


0.1155 


235 



The quality of the NNLO fit to the major data sets is shown in Table 2, together with that 
of the NLO and LO analyses. For each of the smaller data sets, e.g. the CDF VT— asymmetry 



27| and E866 Drell-Yan asymmetry pO], the x pei' degree of freedom is about 1 per point. 



For all the DIS data sets the numbers are quoted for statistical and systematic errors added 
in quadrature.^ The quality of the fits to the individual data sets is satisfactory. For the 
E605 data the systematic errors are quoted in a slightly ambiguous manner, and are generally 
subdominant, and so we fit to statistical errors alone. Hence the quite large fhis case. 

Our treatment of the correlated systematic errors for the Tevatron jet data has been discussed 
in Ref. In Fig. 1 we also show the LO, NLO and NNLO descriptions of the F2 data in a 
few representative x bins. 

For the fit to the DIS data the story is relatively straightforward. In general, as we proceed 
from the LO^NLO^NNLO analysis, there is a sequential improvement in the quality of the 
fit (with a couple of dissenting data sets), with the incremental improvement being smaller at 
the second stage. The LO fit fails because the evolution is a little too slow at small x and also 
too slow, and not quite the correct shape, at large x. This is because there are higher-order 
but small-x-enhanced terms missing in the quark-gluon splitting function and higher-order 
but large- x-enhanced terms missing in the coefficient functions. This can not be completely 
countered even with a very large value of the coupling constant and a large small x gluon. 



discussion of the influence of the correlated systematic errors for the HERA data may be found in the 
Appendix of M. 
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Table 2: Quality of the NNLO fits to different data sets, together with the corresponding 
values for the NLO and LO fits to the same data. 



Data set 


No. of 

data pts 


LO 


ATT / \ 

NLO 


ATATT / \ 

NNLO 


ATATT / \ T 

NNLOJ 


HI 


400 


493 


382 


385 


381 


ZEUS 


272 


296 


254 


264 


275 


BCDMS yup 


179 


179 


193 


173 


175 


BCDMS jjd 


155 


206 


218 


231 


226 


NMC fip 


126 


131 


134 


121 


120 


NMC /id 


126 


98 


100 


88 


86 


SLAC ep 


53 


83 


66 


59 


65 


SLAC ea 


54 


no 

98 


56 


59 


61 


E665 fip 


DO 


oU 


Oi 


00 


oy 


E665 


53 


62 


61 


62 


64 


CCFR F^^ 


74 


116 


85 


80 


78 


CCFR Fg^^ 


105 


126 


107 


110 


126 


NMC n/p 


156 


157 


155 


150 


148 


E605 DY 


136 


282 


232 


220 


268 


Tevatron Jets 


113 


123 


170 


186 


124 


Total 


2097 


2500 


2328 


2323 


2341 



The improvement in the NNLO fit is due to the inclusion of these small and large x terms at 
NNLO, and this also allows a smaller coupling constant than for the NLO fit. However, the 



improvement in the quality of the fit to DIS data is not as dramatic as that seen in Ref. ||15 
This is partially due to our raised W"^ cut - NNLO working a little better in the now excluded 
10 GeV^ < < 12.5 GeV^ range. However, the improvement at small x is not as great, 
and this seems to be due to the fact that the fit actually preferred the previous estimate of 
the NNLO splitting functions which had quicker small x evolution. This may be a sign that 
yet higher orders are important in this regime. For the Drell-Yan data the pattern is the same 
as for DIS data, the improvement at each successive order being correlated largely with the 
reduction in the coupling constant, since as we discovered in the Drell-Yan data prefer a 
lower value of a5(M|). 

The situation for the Tevatron high-i^y jet data is rather different and very interesting. As 
we see from Table 2 there is actually a distinct deterioration in the quality of fit to the jet 
data as we go to higher orders, and at LO we obtain almost as good a fit as seems possible 
given the scatter of points. This is a very surprising result (we are not aware it has been 
noticed previously), but actually quite easy to understand. The way in which the different 
parton scattering combinations usually contribute to the jet rate is shown in Fig. 2 of |31]. 
At the high Et end of the jet measurements the rate is mainly determined by quarks and is 
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given roughly by a'^iqq + qq)- In a standard fit this contribution is usually determined very 
accurately and the normal excess in data must be accounted for by an increase in the g{q + q) 
contribution, i.e. in the high x gluon. Hence the better description obtained in the special 
CTEQHJ and MRST2001J ^ sets. However, at high x the quark coefficient functions are 
enhanced at large x at higher orders, behaving roughly like a5(ln(l — — x))^ at NLO, and 
like a|((ln(l — x)Y)/{l — x))_^_ at NNLO. These positive contributions lead to a smaller quark 
distribution at high x, as one goes to higher order, in order to fit the DIS data. Also, as we 
have seen, the coupling required by the best fit decreases with increasing order. The net effect 
is that o.%{qq + qq) at high x is much (up to 40%) larger at LO than at NLO. At x ~ 0.1, which 
is relevant for the lower Et end of the jet data, where the rate is dominated by the g{q + q) and 
gg contributions, the gluons are a little smaller at LO due to the need to have more gluons at 
very small x, and this compensates for the increased as- Overall, this results in the shape of 
the parton distributions at LO being very suitable for the description of the Tevatron jet data 
without any strong constraints (the gluon at LO is naturally not too small at high x). The 
comparison of the LO and NLO descriptions of the CDF jet data is shown in Fig. 2, and as 
one can see the shape and normalization of the LO theory prediction compares very well to the 
data without recourse to the large shifts in data from correlated systematic errors. As shown 
in Table 2, this translates into a huge reduction in x^- We find this result of a very good match 
at LO, and a big change in going to NLO, very intriguing. This observation may possibly be 
relevant for analyses which extract the jet rates using LO Monte Carlos with LO partons.Q 

Conversely at NNLO the quarks are slightly smaller at high x than at NLO and this means 
that there must be slightly more gluon needed at high x for a good fit to the jet data. Within 
the context of the global fit this means that at NNLO the fit to Tevatron jet data is slightly 
worse and there is an even stronger constraint on the high x gluon than at NLO.[] This increase 
in the fo^' the jet data at NNLO means that the overall improvement in the quality of the 
global fit compared to NLO is only marginal. 

As in our NLO analysis P], we have investigated the sensitivity to as{M^). The results 
are similar, with an approximate error of 0.002 again being attributable due to the quality of 
the fit. The profile is similar to that for the global fit in Fig. 16 of ||^, except that the 
minimum is now at q;5(M|) = 0.1155 and the increase in is slightly steeper for increasing 
as{M^) rather than decreasing as{M^). Indeed the profiles are rather similar to those in |^ 
for each data set. As in where more details of the procedure can be found, we have also 
obtained a special set of partons (denoted NNLOJ) by forcing a very good fit to the jet data. 
The quality of this fit is shown in Table 2. As at NLO the coupling is required to rise to obtain 
the best (NNLOJ) fit, but in this case only to a5(M|) = 0.118. The shape of the NNLOJ 

^We have checked that the good jet fit is a feature of LO partons in general, not just this particular set. 

^In jl^ we did not incorporate independent fits to the jet data at each order, i.e. the same high x gluon 
constraint was used for the NNLO as for the NLO analysis and that at LO was different only by a constant 
K-factor. This is responsible in part for some of the differences between the results in []l5| and this paper, e.g. 
the higher value of as{M'^) used at LO and perhaps the smaller improvement in the quality of the fit at NNLO 
compared to NLO. 
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high-x gluon is very similar to that at NLO, seen in Fig. 15 of |p, and the way in which the 
global fit is affected is similar to NLO, i.e. the main casualty of the improved jet fit is the 
Drell-Yan comparison. 

We now consider the impact on the parton distributions. In Fig. 3 we show the gluon 
distributions for each of the analyses at several values of . At the starting scale of = 
1 GeV^ both NLO and NNLO gluons are negative for very low x with the NLO lying below 
the NNLO. As increases, the NNLO gluon evolves relatively more slowly so that by = 
5 GeV^ the situation is reversed, both gluons being positive with the NLO gluon now above 
the NNLO one. Also shown are the NNLO gluons from the fits using the extreme allowed 
forms of the splitting functions. We note that the spread between these gluons is somewhat 
reduced, particularly above x = 10^^, compared to the range obtained in Ref. |]T3[ as a result 
of the improved reliability of the NNLO splitting functions of Ref. [0. We also note that the 
NNLO gluon is generally closer to the NLO gluon at very low x compared to the case in |jT5[ 



This is due to the decrease of the NNLO quark-gluon splitting function at low x compared to 
the previous expectation, therefore allowing more small x gluon at NNLO. As in Ref. |15] the 
negative feature of the gluon at low x and low is not, in itself, a matter for concern since it 
is not a physical observable, being a scheme dependent quantity. 

In Fig. 4 we compare the parton distributions found in the NNLO fit to those obtained in 
the NLO analysis p. We plot the NNLO/NLO ratios for the gluon, and the up and down 
quark distributions, at = 10 and 10^ GeV^. In going from the NLO to the NNLO analysis, 
several features of Fig. 4 are worth noting. The decrease of the quark distribution at high x 
and the increase at low x (at the lower value) reflect the behaviour of the NNLO quark 
and gluon coefficient functions respectively. We have already seen above that for > 5 GeV^ 
the NLO gluon has overtaken the NNLO gluon and this situation remains as the evolution 
continues. The rise at large x for the NNLO/NLO gluon ratio is a consequence of the need for 
more gluon at high x to fit the Tevatron jet data and also of momentum conservation - the 
positive NNLO quark-gluon splitting function at low x means less gluon is required here and 
more is available elsewhere. This enhancement of the gluon for x above about 0.1 leads to a 
particularly improved fit to the NMC DIS data at NNLO. With the evolution to high Q^, the 
NNLO effects play less of a role, and the NNLO coefficient functions have very little impact. 
However, the relatively smaller low-x NNLO gluon introduced at low produces a slower 
evolution of all partons at low x, which can be seen at = 10'^ GeV^. 

Computing NNLO predictions for the longitudinal structure function Fl requires a knowl- 
edge of the 0{ag) coefficient functions. In Ref. ]15|] we used the partial information of the first 
four even moments together with the x — > behaviour to estimate analogous expressions to 
the 0{ag) splitting functions of Ref. jl^. Here we continue with that approach to compare the 
NNLO prediction for Fl with the NLO one. In fact the two additional moments lead to the 
most likely coefficient functions being essentially unchanged from those used in ||15[ (a lucky 
accident), but we note that the band of uncertainty on these coefficient functions is much larger 
than on the splitting functions. Truly predictive results await a more complete analysis. In our 
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recent NLO analysis we drew attention to the unphysical behaviour of the NLO prediction 
for Fl at low and very low x, shown in Fig. 5, suggesting that the negative feature may 
be a hint that NLO cannot provide the entire description in that region. Now we see that 
the NNLO prediction is indeed more satisfactory, with the physical observable Fl remaining 
positive at = 2 GeV^ down to x = 10~^. Even this qualitative conclusion is not definite 
with the present level of uncertainty. More reliable is the fact that the NNLO prediction is 
below that at NLO at high Q^, since this relies less on the coefficient function and more on the 
gluon evolution. 

The improvement in the description of F2 going from NLO to NNLO is not as dramatic as 
that going from LO to NLO. This could be interpreted as evidence for some convergence as we 
go to higher orders. One phenomenological question that has been raised is whether the data 
which lie below our cut of W'^ = 12.5 GeV^, which is imposed to remove possible higher twist 
contamination, could be consistently described by NNLO effects. At large x the data below this 
cut show a marked rise above the NLO description - as much as 40 — 50% at low - suggesting 
a sizeable higher twist contribution. Fig. 6 shows that the differences between NNLO and NLO 
at large x are indeed significant, and are always in the correct direction. However, they are not 
at the level needed to account for the magnitude of the observed discrepancies. 

We have also investigated the predictions for the W and Z production cross sections at the 



Tevatron and the LHC, and compared them with those obtained previously in Ref. fl^. For 



purposes of comparison we have kept all electroweak parameters, branching ratios etc. at the 



same values as in [p3| , piql . The results are shown in Table 3. As one can see the difference 
between the 2001 and 2000 predictions is always of the order of 0.5% or less, showing that 
neither the new data in the fit nor the changes in the NNLO splitting functions have led to 
any significant changes.^ However, the error band on the prediction due to the uncertainty on 
the splitting functions has been significantly reduced due to the better estimate of the NNLO 
splitting functions. At Tevatron energies the difference between the 'slow' and 'fast' NNLO 
evolution predictions is at the per mille level, while at the LHC the difference is still less than 
1%. 

In summary, we have have carried out NNLO and LO global analyses of data which include 
the most recent data on DIS from HERA together with constraints from hard scattering data. 
Comparisons were made with our recent NLO analysis for the the structure functions and 
for the partons themselves. There is improvement in the general quality of the description of 
the data as one goes to higher orders, though the improvement in the global fit at NNLO is 
not as great as one might expect, and there is also the surprising result that the LO partons 
give an excellent fit to the Tevatron high-ET jet data. The partons at NNLO change by up to 
15% for > 10 GeV^ compared to those at NLO, with the largest change being in the gluon, 
which decreases for x less than about 0.05 but increases for x > 0.1. The quarks increase for 
X < 0.01 at low due to the negative gluon coefficient function, but decrease at lower x for 



^The predictions in Q were made using the 'hybrid' method of combining NLO partons with the known 
NNLO coefficient functions. 
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Table 3: Predictions in nb for W and Z production cross sections at the Tevatron and LHC, 
compared with those of MRST2000 P . 









MRST2000 


MRST2001 


NLO 


Tevatron 


■ 


aw 


2.39 


2.41 




Be+e- 


■ (Tz 


0.219 


0.220 


LHC 


Biu ■ 


aw 


20.5 


20.6 




Bi+e- 


■ (Tz 


1.88 


1.89 


NNLO 


Tevatron 


Biu ■ 


(^w 


2.51 


2.50 




Bi+i- 


■ 


0.230 


0.230 


LHC 


B^y ■ 


aw 


19.9 


20.0 




Bi+i~ 


■ 


1.84 


1.85 



higher because of the slower evolution due to the smaller gluon. They also decrease at very 
high X due to the positive effect of the NNLO coefficient function. 

We note that these are the first generally available sets of partons at NNLO and should add 
impetus to the large activity presently under way to compute NNLO corrections to a variety of 
hard processes. In particular, the computation of the 'missing' NNLO coefficient functions for 
the differential Drell-Yan and large- i^r jet cross sections used in the global fit to constrain the 
large-x sea and gluon distributions respectively is now a matter of some urgency. The NNLO 
partons discussed ini this paper can be found at |http : / /durpdg . dur . ac . uk/hepdata/mrs| 
together with the recent NLO partons (MS scheme) from a similar analysis [0. For completeness 
we have also added the NLO partons for MRST2001 in the DIS scheme and the LO partons 
discussed in this paper. 
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Appendix 

We use a very similar charm prescription to that adopted in [|l^, which is a generalization of 
the Thorne-Roberts variable-flavour-number scheme (VFNS) |Q to NNLO. As explained in 
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||15|| , this has to be a model, since the NNLO fixed-flavour-number scheme (FFNS) coefficients 
are not known. In principle we should be more sophisticated and (i) use a correct matching 
procedure for the NNLO VFNS coefficient function (possible at present), (ii) take into account 
the small discontinuity of the charm parton distribution at NNLO (also possible at present), 
and (iii) also take in account the effect of the NNLO FFNS coefficient function (not possible at 
present). Since the full NNLO prescription is missing, and since the NNLO splitting functions 
are approximate themselves at present, we adopt a simple procedure at NNLO. In JlSl we 
just used the massless C(a|) gluon coefficient function multiplied by the velocity of the heavy 
quark /?=(! — Am'jj z / {Q'^ {1 — z))^'^, which has the correct threshold behaviour but ^ 1 at 
high Q^, as the NNLO coefficient function at all Q^. Here mn is the mass of the heavy quark. 
This has the shortcoming that it leads to a (small) negative contribution at low and small 
X since the the massless C(a|) gluon coefficient function is negative at small x, whereas we 
would expect the NNLO FFNS contribution to actually be positive in this region. Indeed the 
low X, low charm data are underestimated using our procedure. Hence in this paper we 
adopt the minimal simple modification and weight by rather than (3, so that this negative 
contribution is minimized. This improves both the fit to charm data and the global fit. This 
simple approach seems to us to be sufficient at the present level of theoretical accuracy. Once 
the splitting functions are known exactly at NNLO we should improve our NNLO charm model. 
In principle one can even produce an approximate NNLO FFNS coefficient function from the 
known small x limit |^ and known threshold logarithms . 
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Figure 1: The description of data for the F2 structure function at a few representative values 
of X obtained in the LO, NLO and NNLO global analyses. 
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Figure 2: The comparison of the quahty of the LO MRST and the NLO MRST fit to the CDFIB 
high-i?T jet data [Q. The open points are the NLO values before the correlated systematic errors 
have been considered while the solid points are those for the LO fit. 
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Figure 3: Gluon distributions obtained from the LO, NLO and NNLO analyses at various values 
of . The three NNLO gluons result from the slow and fast extremes of the splitting functions 
together with the average of these. 
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Comparison of MRST2001 NNLO and NLO Partons 




Figure 4: A comparison of the gluon, up and down distributions obtained in the NNLO analysis 
with those obtained in the NLO fit 0]. The comparison is shown for two values of Q^. 
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